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I. INTRODUCTION 

The problem of pattern formation is a multi-disciplinary challenge of great technological and scientific interest, and 
has been studied extensively in recent years by physicists, chemists, biologists, materials scientists and others 
Much of the work to date has been been devoted largely to general aspects of pattern- forming instabilities, pattern 
selection, waves and fronts, and direct numerical simulation of the idealized deterministic equations modelling the 
phenomena under study Q. In this approach to modelling, it is reasonable to assume a continuum and coarse-grained 
description for the dynamic variables (degrees of freedom) and the model equations are frequently, though not always, 
of the reaction-diffusion type. Thus, one begins by writing down some specific deterministic model equations. This 
may be regarded as the phenomenological approach to the problem. However, a first principles derivation of the 
dynamical equations is obtained by removing the fast or short wavelength degrees of freedom from the microscopic 
description, and this process leads to unavoidable noise terms, representing the effects of internal fluctuations, which 
result from this small-scale elimination step [3|. Additionally, in the the study of complex phenomena, there is 
typically no precise knowledge of many of the microscopic details, nor of the initial or boundary conditions needed 
to provide a complete description of the problem. The dynamics may also evolve in a medium (such as a background 
fluid) which provides external perturbations, environmental noise and unpredictable disturbances. Nevertheless, one 
is interested in an explicit understanding of the system at long wavelengths. For all these reasons, it is natural to 
consider the influence of random noise on otherwise deterministic models and to study stochastic reaction-diffusion 
equations. 

One of the simplest models of biochemical relevance leading to spatial and temporal patterns when diffusion is 
included is that due to Gray and Scott [J] . Numerical simulations of deterministic systems have revealed a surprisingly 
large set of hitherto unknown complex and irregular patterns [Bj]. Because of the above considerations, and in regards 
to these recent findings, three questions that immediately come to mind are (a) how do the fluctuations affect the 
stability of an established pattern? and (b) what are the emergent properties, due to fluctuation effects, of such a 
system at long wavelengths? and (c) how do the deterministic and stochastic effects compete? The purpose of this 
paper is to formulate carefully an analytic answer to the second question posed in (b). The latter will be carried out 
with the help of the dynamic renormalization group (RG). The questions raised in (a) and (c), of how noise influences 
pattern selection will be investigated numerically and will be presented elsewhere. Striking numerical evidence for 
noise controlled pattern self-replication is discussed in [||, where the replication rate is maximal for an optimal but 
small noise intensity. 

We therefore consider a stochastic version of the Gray-Scott model [4j, defined by the following system of stochastic 
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partial differential equations: 

^-V = XUV 2 - fiV + D V V 2 V + ^u(x, t) 

= u - XUV 2 ~vU + D U V 2 U + J7„(x, t). (1) 

ot 

In the absence of noise, (JTJ) coincides with the Gray-Scott model which is a variant of the autocatalytic Selkov 
model of glycolysis, corresponding to the following chemical reactions: 
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V (x, i) and (x, t) represent the concentrations of the chemical species U and V, and are functions of d-dimensional 
space x and time t. A is the reaction rate, P and Q are inert products, fi is the decay rate of V and v is the decay 
rate of U and u a is the constant feed rate. In the reaction rate parameter was simply set equal to unity: A = 1. 
Here we retain it as a free parameter, as we will see that it undergoes a nontrivial renormalization. A non-equilibrium 
constraint is represented by a feed term for U. The rate at which U is supplied is positive if the concentration of U 
drops below an equilibrium value and negative if it exceeds it. The equilibrium U concentration is uq/v, where u$ is 
the feed rate constant. The chemical species U and V can diffuse with independent diffusion constants D u and D v . 
All the model parameters are positive. 

Any real chemical system is subjected to random fluctuations. We can include such effects in the GS model by 
means of noise terms, that can in principle be additive or multiplicative. In this work, we have chosen to investigate 
the influence of additive noise alone as a initial approach to the full stochastic problem in which both types of 
fluctuations can be simultaneously present. This is intended therefore as a first step towards incorporating fluctuation 
effects. This choice is influenced in part by the technical aspects of the dynamical renormalization group. From 
this perspective, additive noise, which enters the dynamics linearly, is technically easier to treat analytically in RG 
calculations. Accordingly, here we have opted to use independent additive colored noises, ?7 t) (x, t) and rj u (x,t). First, 
spatially-correlated noise, and then temporally-correlated noise are considered. A few cautionary remarks regarding 
our incorporation of noise are in order. We emphasize that here we are adopting a simple phenomenological approach, 
wherein noise is added to the idealized deterministic equations. But, stochastic equations can also be derived from 
first principles. We refer to the methods that take a classical master equation, purporting to encode the precise 
microchemistry, to a continuum field theory @, H, • This " second-quantized" formalism requires detailed knowledge 
of the microscopic master equation, and the end result of this method yields the stochastic PDE for the coarse- 
grained degrees of freedom and the (internal) noise correlations. A technical discrepancy comes up in applying 
phenomenological reaction-diffusion equations with real additive noise to pair- reaction kinetics (i.e. V + V — > 0), 
since in this particular case of particle annihilation, the SPDE derived from the fundamental microscopic master 
equation is complex and contains imaginary noise [l0|. However, for standard Gribov processes (particle clustering 
reactions), the first-principles method docs lead to a real SPDE with real additive noise. Thus, the ad- hoc method is 
adequate for handling these cases. For external or environmental noise, the phenomenological strategy adopted here 
is also adequate since external fluctuations are typically specified only at the coarse-grained level. 

The remainder of the paper is organized as follows. In the next section, we specify the noise properties directly 
in Fourier space and derive the scaling laws that the model parameters in (J]} obey when the stochastic equations 
are scale invariant (as they are in the neighborhood of any RG fixed point) and the general scaling form of the 
correlation functions of the composition fields is deduced in Subsection III Al This information will be needed when 
we apply the renormalization group in Section [Ml The fixed points and associated exponents are solved for and the 
RG flow is represented in a reduced two dimensional parameter space. The global flow patterns clearly illustrate the 
d-dependence of the flow, the role of the associated critical dimension, as well as the phenomena of crossover. In 
Section [TV] we discuss the emergence of effective three-body reaction terms not present in the original GS model that 
arise as a direct consequence of fluctuations and non-linearity. Modifications and changes in the RG results due to 
the presence of temporally correlated noise are very briefly discussed in Section |Vj Conclusions and discussion are 
presented in Section PVTl 
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II. NAIVE SCALING PROPERTIES 



We can deduce the naive scaling laws that must hold if the stochastic equations ([I]) are to be form-invariant under 
a basic re-scaling of both space and time. To do so, we specify the properties of the noise. In Fourier space, the noise 
correlations are given by 

(r) v (k,uj)ri v (k',u')) = 2(2ir) d+1 A v k~ y *6 d (k + k')5(uj + uj') 
(t ?u (k,w)»?u(k',w / )) = 2(2Tr) d+1 A u k-y-S d {k + k')S{u; W) 

(r hl (k,u;)r )u (k',uj')) = 0. (3) 

We assume all cumulants except the above to vanish, thus the noises are individually Gaussian, centered around 
zero and mutually uncorrelated. The independent noise amplitudes A v , A u > are positive definite and the noise 
exponents y v and y u are real free parameters. 

Space, time and the concentrations [/, V , are re-scaled by the following transformations: 

x — > s x 

t -> s~ z t 

V -> s~ x "V 

U s~ x ^U. (4) 

Here, s > 1 is a convenient scale factor and z is the so-called dynamic exponent. We initially allow for each 
concentration field to respond independently under the re-scaling; this is reflected by the two a-priori independent 
roughness exponents \ v and Xui respectively. However, the fact that the reaction rate A appears simultaneously in 
both the U and V equation requires that Xu = Xv = X: as is easy to check. There is just a single roughness exponent 
to deal with. 

Subjecting the equations (fT]) to this transformation, we find the model parameters scale naively (classically) in the 
following way, 
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A. Correlation functions 



Both the dynamic and roughness exponents govern the scaling form of the correlation functions of the concentration 
fields. To see this, consider the effect of a scale transformation on the y-field which we write more explicitly as 

V{sx,s z t) = s x V(x,t), (6) 

which holds in the scaling regime (when the system is near one of its RG fixed points). Then the correlation function 
will scale as 

(V(x,t)V(0,0)) = »- 2 *(I / (»x,»"t)l'(0,0)) 
where the scaling function $ obeys the following asymptotic limits (TlT | 



lim $(w) = const. lim = const. xu 2x ^ z . (8) 
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Thus, knowledge of both scaling exponents [z, x) is required in order to predict the long-wavelength and long-time 
correlations of the fields. Similar considerations apply to the (UU) and cross-correlation (VU). In fact, the scaling 
of these latter two correlations will obey a relation similar to (O, except for a possibly different scaling function. 
Nevertheless, the limits in ([5]) are valid for all scaling functions, although the constants appearing there can be 
distinct. The tool best suited for the calculation of the required exponents is provided by the renormalization group, 
to which we next turn. 



III. RENORMALIZATION GROUP (RG) ANALYSIS 



Since we are interested in the emergent scaling properties of the model, we focus on the hydrodynamic limit, that 
is, the long-time t — > oo and long-distance x — > oo limits. The dynamical RG [HI, [l3| is a powerful tool for computing 
and analyzing the asymptotic properties of out-of-equilibrium stochastic systems. 

In order to calculate the scaling exponents in the asymptotic limit, the linear and non- linear parts of the equations 
of motion can be re-organized so as to allow a perturbative calculation in powers of the non-linear parameter A. After 
Fourier transforming, the coupled stochastic equations of motion (TTJ) can be rewritten as 



2tt 



V(k,u>) = G v0 (k,Lu)ri v (\s.,u) 

cZ d ki duji d d k 2 duj 2 

G„ (k, w)A 

U(k,cu 

G„ (k, w)A 



(27r) d 2tt ( 27r ) 



V(k x , u3i)V(ka, w 2 )t/(k - ki - k 2 ,u> — wi - uj 2 ) 



G u0 (k, uj) Vu (k, u) + G uQ (k, Lu)(2ir) d+1 6 d (k)S(uj)u 
d d ki dioi d d k 2 dw 2 



(9) 



V(k 1 ,uji)V(k2,Lu 2 )U(k - ki - k 2 , w - lo\ - w 2 ) 



(27r) d 2tt (2n) d 2vr 

where the bare propagators (or, response functions) G v0 and G u0 are defined by 

G„ (k,u) = 



H - 



D v k 2 — iuj 
1 



v + D u k 2 — iuj 

The wave-number modulus is denoted by k — \k\. Define the effective propagators by 

V(k,u) = (3f tl (k,u)» 7o (k,w) 

U(k,u) ee G u (k^){ Vu (k,u;) + (27r) d+1 S d (k)S(io)u a }. 



(10) 



(11) 
(12) 



Substituting these into the set of integral equations ([9]) yields a set of equations that can be solved iteratively to any 
order in A. To this end, it is best to handle the expansion via diagrams, and the above equations © are represented 
diagrammatically as follows: 



Vv 



Vv 




The directed double straight (wiggly) line symbol represents the effective ^-propagator ({/-propagator) defined above: 
straight and wiggly lines correspond to the V-sector, and [/-sector, respectively. The single directed lines stand for 
the bare propagators (|10|) . It is understood that each directed line carries wave number k and frequency to, which 
is conserved at all vertices, as can be seen from We do not write down this dependence explicitly on the graphs 
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in order to avoid clutter. The non-linear coupling terms, or vertices, are represented by the directed four-pronged 
symbols, of which there are two, and these are denoted by the encircled numerals ® and (2), respectively. The 
zero-mode uo in the second equation in (J9j) is indicated above by a (^). The perturbation expansion can now be 
performed efficiently without having to carry along lengthy and tedious algebraic expressions. For further details on 
using graphical methods to solve stochastic differential equations, see the Appendices of [141 ] . The calculation of the 
propagator then follows from a graphical iteration of these expressions, amputating one noise factor from the legs, 
followed by an averaging over the remaining noise factors. The results at one-loop level are indicated as follows: 




+ 2 loops 



Note that only the U propagator receives corrections at one-loop, whereas the V propagator is unchanged to this 
order. Furthermore, the v- noise is what "drives" the renormalization of the [/-propagator. The loop expansion for 
the vertices is indicated below. One takes the vertices and iterates, to a desired order in the coupling A, by replacing 
directed double lines using the above equations of motion. Then, three noise factors (two ry^'s and one r] u ) are 
amputated and the resultant graphs are averaged over the remaining noises. A combinatorial factor of 4 appears after 
counting all possible noise contractions. We thus obtain, 




+ 2 loops 



+ 2 loops 



The noise spectral functions receive no corrections at one loop order. We point out that although A is employed 
as a formal expansion parameter, the bona-fide and dimensionless perturbation parameter g (see (|14[) below) is a 
combination of this, a noise amplitude, a decay constant, a diffusion constant and a short distance (or ultraviolet) 
cutoff needed for convergence of the integrals. Moreover, we expand in loops rather than in powers of g. Physically, 
the loop expansion is an expansion in powers of the noise amplitude [la ] . We note that at leading order in loops, the 
perturbation expansion indicates that only the U field propagator and the non-linear coupling A receive corrections. 
The noise amplitudes and the V propagator are unchanged at one-loop order. At two-loop order, all propagators, the 
noise amplitudes as well as the vertices, do receive corrections. However, there are some important issues at leading 
order that need to be understood. 

The dynamical RG transformation is carried out in two steps (l3| : 

1. High momenta components, in the momentum shell Ae~ Al < |k| < A, at and below the cutoff A, are integrated 
out. Note: s = e l . 



2. A change of scale restores the cutoff to the value A. 
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Further details of the calculation are given in Appendices A and B. Working with infinitesimal parameter SI gives 
rise to differential equations governing the RG flow in parameter space: 
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(13) 



Kd = Sd/(2ir) d where Sd is the surface area of a ci-dimensional sphere. As already noted in the corresponding graphs, 
at one-loop order, only two out of the total of eight model parameters run with scale, namely, the decay rate v, of 
the U field and the nonlinear coupling A. Note moreover that at this order, the corrections are driven by the w-noise, 
but not the it-noise. Despite the relative simplicity of these equations, non-trivial RG fixed points and flow already 
result at this leading order perturbation. 



A. Fixed Points and Dynamic Scaling 

We determine the fixed points and associated exponents implied by the above equations (TT5|) . We can analyze the 
RG flow in a reduced two-dimensional parameter space by introducing the pair of dimensionless couplings defined by 

a d-v v -2 

9 = XA v K d =— , (14) 

uD v 

A 2 

h = (D v + D u ) — , (15) 

for then the RG equations (TIB)) can be written as 

v = u{z + g), (16) 

we do not write out the remaining trivial RG equations. Note we have taken fi — from the outset. This choice is 
necessary in order to obtain a non-trivial g ^ fixed point solution of (|13p . Here the over-dot stands for the derivative 
dj dl. In terms of g and h, the nontrivial part of the RG flow is governed by the pair of equations 

9 = g (yv-d + 2-g--^-) , (18) 



l + h / 

h = -(2 + g)h. (19) 

The fixed points can be solved for by looking for all zeroes of the pair (|18I19[) . to be denoted as (g* , h*) while the 
associated fixed point exponents are obtained by substituting the solutions of g = and h — into (|16I17P ; solving 
for the zeroes of this latter pair then yields the exponents, denoted by (z*,%*). 

We first search for all non-trivial fixed points. These correspond to g* ^0, since g is proportional to the non-linear 
coupling A. There are two cases to be distinguished. 



7 



Case(a): g* = —2 and 1 + h* = 8/(d — y v — 4). The associated exponents are z* = 2 and x* — (2 — d + y v )/2. 
Note that the combination (z* — 2\* — d + y v ) = vanishes identically. Thus, if we choose uq — 0, and y u = y v , 
all remaining RG equations (fl~3|) are stationary. Note for this fixed point solution, the stochastic GS model is in the 
same universality class as the linear Edwards- Wilkinson model [14| . 

Case(b): h* = and g* — (y v — d + 2)/5. The exponents associated to this fixed point are z* — —(y v — d + 2)/5 and 
X* = (Uv — d + 2)/2. Again we choose uq — 0, then all remaining RG equations are automatically stationary except 
that the noise amplitudes decay to zero on approaching this fixed point: A v (l) = A v e l ( z ~ 2x ~ d+Vv ) — ► as I — ► oo and 
similarly for A u (l). This limit holds provided that d < 12 + y v , and d < 12 + y u for the limit of A u . 

Since the model (JTJ) is defined as having all non-negative parameters, we must exclude the solution in case(a) as 
being physically spurious. This leaves us with the nontrivial fixed point of case(b), which is an attractive fixed point 
and is indicated by the symbol A in the flow graph; refer to Fig. (jT]) . 

There is one trivial fixed point: this is a saddle point, as indicated by S in the flow graph Fig.fl]). This corresponds 
to g* = and h* — 0. The exponents are found to be z* = and x* — [Vv — d)/2. This is a fixed point solution to 
all equations in (fT5)l provided y v — y u and again upon setting uq = 0. 

B. RG flow 

The RG flow in the regions of parameter space surrounding the fixed points is obtained by numerically integrating 
the differential equations (|18I19|) for various choices of distinct initial conditions chosen from within the basins of 
attraction or repulsion of the fixed points. The results are shown in the Figs. (|1I2I3|) . The topology, the direction 
of the flow and the fixed point stability are controlled by the single parameter epsilon e = 2 — d + y v . In deriving 
cases (a) and (b) above, we tacitly assumed that e > 0. The corresponding flow is plotted in Fig.fl]). The origin is a 
saddle point and there is one asymptotically stable fixed point, as shown there. For very small initial g, the effective 
dynamics will flow towards S and then be repelled to A. Thus, the system exhibits crossover. In the vicinity of S, 
the correlations of the concentrations therefore scale as (r = |x|) 

{V{^t)V{0,Q))~r^- d \ (20) 

whereas on approaching the point A, they scale according to 

<V(x,t)V(0,0)) ~ r ( ^- d+2 \ for oo, (21) 
- t- 5 , for t -> oo. (22) 

If e = 0, there is no fixed point at nonzero g. The origin (0, 0) changes from a saddle point to an attracting sink 
as e — > 0. The critical dimension is defined by e = 0, and is given by d c = 2 + y v . This is the dimension below 
which fluctuations are relevant. For example, for white noise, d c = 2, which is also the critical dimension of the 
Edwards- Wilkinson model [lj]. 

For e < 0, there is only the trivial fixed point A at the origin. Both these cases are shown in Fig. ^ and Fig. ([3]), 
respectively. In the neighborhood of this point, the correlations again scale as in (|2"D|) . 

IV. INDUCED THREE-BODY REACTIONS 

Consider the graphical representation of the stochastic equations of motion written in ^ and which are displayed 
below Eq. (fl"2"l) . In arriving at the perturbative expansion for the vertex representing the chemical reaction UV 2 in JT]) , 
the steps used in the graphical method involved replacing two of the vertex legs with the exact expression for V(k, oj) 
and the third remaining leg with the exact expression for £7(k, u>) using Eq.Q. This step must be iterated once more 
in order to obtain all one- loop diagrams. At this intermediate stage however we generate a set of tree diagrams (i.e., 
containing no loops) with all legs terminating either in an rj v or in an r\ u noise factor. Now, the effective UV 2 vertex 
results after amputating (i.e., removing) two 77^ 's and one r] u noise factor from each tree diagram and then averaging 
over the remaining noise factors. The averaging step makes use of the noise correlation functions specified in Q. The 
combinatorial factor for each diagram is found by counting all possible noise contractions that give rise to it. In this 
way, we arrived at the one-loop expansion for A depicted diagramatically above in Section [Ml 

However, at the step leading to the set of tree diagrams, one can instead amputate two r; u 's and one rj v noise from 
the tree diagrams and then average over the remaining noises or, alternatively, amputate three rj v factors and then 
average. The sequence of noise factor amputations can lead to different allowed reaction kinetics. Thus, at one-loop 
level, in addition to the Gray-Scott UV 2 reaction, one also encounters the following (induced) effective reactions 
involving a triple product of chemical concentrations: 
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0.25 



FIG. 1: A projection of the RG flow in Eqs. (|18I19|I . Two fixed points: A is attractive, S is a saddle point. This is plotted for 
e > 0. 




0.25 



FIG. 2: A projection of the RG flow in Eqs. (|18I19[1 . One fixed point: A is attractive. This is plotted for e = 0. 



U 2 V 



V 6 



(23) 



It is important to realize that these are alternative but exclusive reactions, in that noise-factor amputation can 
produce one reaction, say (U 2 V), or the other (V 3 ), but not both simultaneously. This is simply because each 
alternative reaction derives from the same tree graph, the only difference is due to the sequence of noise amputations 
and subsequent averaging steps. Note it is not possible to amputate three rj u factors at one loop level (this would have 
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FIG. 3: A projection of the RG flow in Eqs. (|18I19|I . One fixed point: A is attractive. This is plotted for e < 

led to a U 3 reaction, which is therefore ruled out from the effective induced dynamics at this order). Both reactions 
in (f2"5]l are proportional to A 2 at one-loop order. Moreover, since at one loop, both concentration fields U and V scale 
with the same exponent x, these effective three-body reactions scale as s 3x and will be relevant (or, irrelevant) if and 
only if the original GS reaction is. Most importantly, the new reactions do not destroy the scale invariance of (Q]) at 
the fixed points. 

Physically, this means that at one-loop order in the fluctuations, corrections are not only induced in the original 
GS reaction, leading to its renormalization (i.e., last equation in (fT3|) ). but that a new effective reaction is generated. 
This new reaction is not present in the GS model (i.e., does not exist at tree level) but is a direct consequence of the 
combined effect of noise and non-linear terms in ([1} . In this sense, the stochastic version of the GS is not renormalizable 
at long wavelengths, because new relevant reaction terms are thereby induced. The large scale chemistry that emerges 
at the attractive fixed point A is GS with corresponding effective (renormalized) parameters plus one of the alternative 
induced reactions written in or in (j2"5| . This is interesting because it indicates that at large scales, the coarse- 
grained theory ([T]) represents an apparently distinct chemistry. 

In terms of chemical reactions, the new reaction vertices would correspond to production of V and or U molecules 
in one of the alternative pathways indicated below. In the case of the vertex U 2 V, we have various pathways 



Thus, for the case of the the induced U 2 V vertex, and for the specific pathway V + 2U — > 2V + U, the effective 
one-loop reaction dynamics at the attractive fixed point is given by the following deterministic equations: 




(24) 



whereas for the case of the vertex V" 3 , 




(25) 



\*UV 2 + (X*) 2 U 2 V + D V V 2 V 
-\*UV 2 - (X*) 2 U 2 V - v*U + D U V 2 U, 



(26) 
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where A* and v* denote the fixed-point values of A and v, respectively. Recall that only these two parameters 
renormalize nontrivially at one- loop order. Moreover, from the fixed point analysis Case(b) given in Sec ((III A[) we 
recall the solution A is obtained for \i = 0, uq = and both the noise amplitudes decay to zero upon approaching the 
fixed point provided d < 12 + y. For an attracting fixed point, e = 2 — d + y > which requires d < 2 + y, which 
is the stronger inequality. Alternatively, for the case of the one-loop induced three-body reaction V 3 , and for the 
pathway 3V — ► 2V + U, the effective one-loop reaction dynamics at the attractive fixed point is given by the pair of 
deterministic equations: 

X*UV 2 - (X*) 2 V 3 + D V V 2 V 

-X*UV 2 + {X*) 2 V 3 - v*U + D u V 2 U. (27) 

The effective reaction dynamics derived above results from the lowest order one-loop perturbative RG calculations. 
If the RG program is carried out to higher order (for example, to two loop order), then these one- loop induced terms 
would have to be taken into account. 



at 
d_ 

dt 



v 



u = 



V. TEMPORALLY CORRELATED NOISE 



Starting from a microscopic description, the elimination of the fast degrees of freedom leading to fl} can in principle 
result in noise with long-range correlations in both space and time. It is therefore of interest to investigate the 
influence of temporally correlated noise on the present phenomenological model. It is straightforward to extend the 
RG analysis to incorporate long range temporal correlations. Since the basic calculational steps are similar to those 
employed above, we will be concise and list only the salient features specific to temporal correlations. 

So, in place of ([3]), we consider Gaussian noise with correlations that behave asymptotically as follows: 

(ri v (k,uj)r} v (k',u>')) = 2(2ir) d+1 A v k-y"Lu- 2e "6 d (k + k')6{L; + Lu'), 

{r) u (k,u>)r} u (k',Lj')) = 2(2ir) d+1 A v k- y "u;- 2e "5 d {k + k')5{uj + u;'), (28) 

where the exponents Vl 6 u control the range of the temporal correlations; the limit of purely spatial correlations is 
recovered by setting 8 V = 6 U = 0. 

Scaling properties of the noise amplitudes are modified as follows: 

A v -> s y«-d-2x+*(i+2e v ) Av 

A . „j/„-<i-2x+2(l+26>„) a 

SX U 7 & -^115 

(29) 

while the remainder of the relations in ([5]) are unchanged. The steps needed to carry out the perturbation expansion 
and associated diagrammatic development of fl]) are the same as before; only now, the noise factors appearing there 
and in the diagrams are those corresponding to (|28p . The calculation of the required loop diagrams and integrals 
follows the same basic steps as outlined in Appendices A and B. Due to (f2"5)) the integration over internal loop frequency 
f2 is much more complicated, though still analytically tractable [IB]. The two non-trivial one- loop differential RG 
equations in (|13p are modified accordingly: 



dv 7T \A v K d A d -V" 

— - -" + csc[(l + 2fl„)-] - ; - A2U+2g 



dl ' LV J 2'(fj, + D V A 

dX 

~dl 



2 X + Z - 4csc[(l - *.,=] W> + D^) { %^ D ^_% + D :Jf 



A, (30) 



and the allowed range of the noise exponent is — \ < 9 V < |. As in the case of pure spatially correlated noise, we can 
analyze the RG flow and fixed points in terms of a convenient choice of dimensionless parameters (taking /i = from 
the outset): 

g = AA^ dCS c[(l + 20„)|] v{J £l 3 y +Mv , (31) 
A 2 

hi = D v — , 
v 

A 2 

h 2 = D u — . (32) 
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This should be contrasted with the pair used in Section IIII At here the temporal correlations do not allow one to 
combine hi and h 2 consistently into a single unified variable (hi + h%) — * h, as we did above. 
In terms of these variables, the RG equations can be written as 

/ ( fe i U+28„ _ 1 \ 

9 = ° [*> - d + 2{l + 2 ^ - 9 - ^ + ^ (i + ^ + mm-^ + m j 

hi = -h x (2 + g) 

h 2 = -h 2 (2 + g). (33) 

The fixed points are solved for by looking for all the zeroes of the triplet of equations in (|33|) . As before, there are 
two non-trivial fixed points and one trivial fixed point. One of the non-trivial fixed points corresponds to g* = —2, 
h^ and h 2 are arbitrary constants, and we discard it for the reasons given earlier (see comments concerning case(a) 
in Section UlI A|) . The other non-trivial point corresponds to g* = (y v — d + 2(1 + 20„))/5 and h\ — h 2 = 0. The 
associated exponents are found to be z* — — (y v — d + 2(1 + 28 v ))/5 and x* = (Vv — d + 2(1 + 29 v ))/2. Previous results 
pertaining to pure spatially correlated noise are immediately recovered by setting 8 V = in these expressions. The 
trivial fixed point corresponds to g* = h\ = h% = 0, with exponents z* = and x* = (Vv — d)/2, results which are 
seen to be identical to those of the purely spatially correlated case. In this case, the flow is governed by the parameter 
e = y v — d + 2(1 + 26 v ). The RG flow is qualitatively the same as before, so in this respect, temporally correlated 
noise leads to no new features at one-loop order. 



VI. CONCLUSIONS AND DISCUSSION 



We have studied the large wavelength and long time limits of the Gray-Scott model subject to random fluctuations. 
We carried this out for additive noise containing long range correlations in space and in time, which leads one 
to consider a set of coupled stochastic partial differential equations. The noise is intended to model in part the 
combination of coarse-grained external fluctuations, environmental noise and also imprecise knowledge of initial and 
boundary conditions. The asymptotic behavior of this system is revealed from applying the dynamical rcnormalization 
group combined with perturbation theory. This behavior is summarized by flows in parameter space, which indicate 
how the parameters of the dynamic model change under coarse-graining; general points in this space represent the 
dynamics effective at the corresponding space and time scale. The fixed points control the effective dynamics at the 
largest scales and corresponding correlation functions (in terms of the chemical concentrations or composition fields) 
are power laws depending on a pair of scaling exponents. It will be noted that the information obtained here is of a 
statistical nature, since here one deals in correlation functions, which are themselves probabilities. These correlations 
provide a measure of the long wavelength properties of the patterns that can form in the presence of fluctuations. 

The RG calculation was carried out to one-loop order. At this lowest order, only two out of eight of the model 
parameters receive corrections, yet there is already a non-trivial infrared stable fixed point and the phenomena of 
cross-over. We have also identified the critical dimension below which fluctuations are relevant. This dimension is 
d c = 2+y v for spatially correlated noise or d c = 2+y v +2(l+20 v ), for noise with both spatial and temporal correlations, 
respectively. In the limit of large scales (after coarse-graining) the properties of the system converge towards a simple 
solution with fractal properties (i.e., exhibiting scale- invariance) as shown in Fig{T]). The existence of this solution 
depends on the sign of e. The non-trivial scale invariant solution is present only for positive values of e, when the 
fluctuations become important and then emergent behavior arises. Moreover, we find that the combined interplay of 
the nonlinearity and the fluctuations lead to effective three-body reaction terms not present in the original GS model. 
At one-loop order in the fluctuations the new reaction vertices correspond to three-body molecular reactions with 
various alternative pathways as expressed in Eqs (pM|) and (|25p . 

The renormalized reaction-diffusion equations therefore can be represented by a chemistry apparently distinct from 
the original GS reaction. For small noise amplitudes, a one-loop calculation should be adequate to capture some of the 
salient features of stochastic reaction-diffusion dynamics. In fact, simple white noise can lead to rather striking effects. 
Preliminary numerical calculations demonstrate clearly that both pattern selection and dynamic pattern replication 
can be controlled by adding white noise to the deterministic GS model These studies raise important questions 
pertaining to the role of noise in both chemical and biological self-organization and the environmental selection of 
emergent properties. 

At one-loop order, temporally correlated noise leads to no noteworthy features in the RG analysis, though that 
conclusion may change at next higher order. Temporally correlated fluctuations can lead to complex scaling exponents 
which signal the presence of hierarchical structures 171 ] ■ Such structures are absent at one loop order. 
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APPENDIX A: PROPAGATOR RENORMALIZATION 



Details of the propagator renormalization are give here. As the diagrams for the expanded propagators show, at one 
loop only the U propagator receives corrections. Reverting back from diagrams to corresponding algebraic quantities, 
we have that 

G„(k,w) = G u0 (k, w ) + A[G tl o(k,^)] 2 / 2A vP -y*G v0 (p,n)G v0 (-p,-n), (Al) 



/ 2A uP -^G v0 (p,n)G v0 (-p,-n), 

Jp,n 



where f ^ = f (£j\d I^°oo 1S an abbreviation for the integration over wave-vector and frequency. The integration 
over wave-vector can of course be further decomposed into an integration over angles and modulus. Eq (|Aip can be 
used to obtain an expansion for the inverse propagator: 

G- 1 (k,w) = G; 1 (k,u;)-A f 2A v p-y*G v0 (p,n)G v o(-p,-n) + 0(\ 2 ). (A2) 

Jp,n 

Referring back to the structure of the bare propagator in (|10|) , in order that the original model JT]) be renormalizable 
in the hydrodynamic limit, the inverse effective propagator must have the form 

G^ 1 (k, oj) = v + D u k 2 — iuj, (A3) 

where v and D u are the effective, or renormalized decay rate and diffusion constants for the U field, respectively. 

Since the one-loop correction integral in (|A2j) does not depend on either u> nor k, we can immediately conclude that 
the diffusion constant is unrenormalized at one loop order, D u = D u while the decay rate renormalization is given by 

v = v + \ [ 2A v p-y»G v o(p,n)G v0 (-p,-n). (A4) 
Jp.n 

The frequency integration over f2 can be performed exactly (e.g., by residues). Finally, eliminating a finite band of 
large wave-numbers below the cutoff A in the Wilsonian fashion then yields 

r A p d-y v -i 

v < =v + \A v K d dp / (A5) 

For an infinitesimally thin wave- number shell (s = 1 + 5, < 5 << 1 ) we pass from an integral to a differential 
relation. After a further re-scaling according to ([5]) we obtain the first differential RG equation as displayed in (fT3"|) . 



APPENDIX B: VERTEX RENORMALIZATION 



Steps similar to those above are involved in the renormalization of the one- loop vertex (or, coupling A). Transcribing 
the vertex diagram in the text back into algebraic quantities, we have that 

A = A-4A 2 f 2A vP -y^G vQ (p,n)G v0 (-p,-n)G uQ (-p,-n), (Bl) 

where we have set all external wave-numbers and frequencies to zero from the outset, in anticipation of the hydrody- 
namic limit. As before, the frequency integration can be performed immediately by means of residues. Eliminating a 
finite band of large wave-numbers yields 

f A „d-v«-i 
A< = A-4A 2 A„A^ / dp P (B2) 

Ja/ s (fJ, + D v p 2 )([i + is + [D u + D v \p z ) 
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For an infinitesimally thin wave- number shell (s = 1 + 6, < 5 « 1 ) we pass from an integral to a differential 
relation. After a further re-scaling according to ([5]) we obtain the last differential RG equation as displayed in Q13p . 
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